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Abstract 


Nomenclature 


The Swann and Pittman semi-empirical 
relationship 1 has been used as a standard in aerospace 
industry to predict the effective thermal conductivity of 
honeycomb core panels. Recent measurements of the 
effective thermal conductivity of an adhesively bonded 
titanium honeycomb core panel using three different 
techniques, two steady-state and one transient radiant 
step heating method, at four laboratories varied 
significantly from each other and from die Swann and 
Pittman predictions. Average differences between die 
measurements and the predictions varied between 17 
and 61 percent in die temperature range of 300 to 500K. 
hi order to determine die correct values of the effective 
diermal conductivity and determine which set of die 
measurements or predictions were most accurate, die 
combined radiation and conduction heat transfer in die 
honeycomb core panel was modeled using a finite 
volume numerical formulation. The transient radiant 
step heating measurements provided the best agreement 
with the numerical results. It was found that a 
modification of the Swann and Pittman semi-empirical 
relationship which incorporated the facesheets and 
adhesive layers in the thermal model provided 
satisfactory results. Finally, a parametric study was 
conducted to investigate the influence of adhesive 
diickness and thermal conductivity on the overall heat 
transfer through the panel. 


A cross sectional area 
F radiation shape factor 
k diermal conductivity 
L height (thickness) 

L' honeycomb core height 
q" heat flux 
r radial coordinate 
T temperature 
z axial coordinate 
e Emissivity 
8 Kronecker delta 

A. height to diameter ratio for honeycomb cell 
c> Stefan -Boltzmann radiation constant 

Subscripts: 
ad adhesive 
f foil 
e effective 
fs facesheet 
g gas 
r radiant 

Acronyms 

ASTM American Society of Testing and Materials 
GHP Guarded Hot Plate 

NIST National Institute of Standards and 

Technology 

SRM Standard Reference Material 

S&P Swann and Pittman 


Aerospace Engineer, Member AIAA 
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Introduction 

An adhesively bonded titanium honeycomb core 
skin panel was considered for the whig structure of die 
High Speed Civil Transport not only for its structural 
performance, but also because of its low effective 
diermal conductivity in order to minimize heating of the 
fuel stored in die wings. Measurements of die effective 
diermal conductivity of typical honeycomb core panels 
using diree different standard techniques at four 
different diermophysical property measurement 
laboratories in die temperature range of 300-500 K 
yielded results diat were significantly different from 
each odier and from die semi-empirical Swann and 
Pittman 1 model for combined radiation and gas and 
solid conduction heat transfer in honeycomb core 
panels. The Swann and Pittman model has been used 
as die standard model for determining heat transfer 
dirough honeycomb core structures duoughout die 
aerospace industry. The significant difference in die 
measured and predicted values of effective diermal 
conductivity and the wide variation among die 
measurements using various techniques required a more 
diorough investigation of heat transfer dirough diis 
specific honeycomb core panel. 

The purpose of the present investigation was to 
determine die effective thermal conductivity of 
honeycomb core panel. This was accomplished by a 
finite volume numerical solution to the combined 
radiation/conduction heat transfer for the specimen. 

The results were then compared with the Swann and 
Pittman predictions and the experimental results to 
determine which set was most accurate. Then, 
modifications to the Swann and Pittman model that 
incorporated the facesheets and adhesive layers in the 
thermal model were investigated and compared with the 
finite volume numerical results. Finally, a parametric 
study was conducted to investigate the influence of 
adhesive thickness and thermal conductivity on the 
overall heat transfer through the panel. 

Background 

Various researchers have investigated heat transfer 
hi honeycomb core structures. Most of the work 
reported in the literature has been limited to panels with 
the honeycomb core brazed to the facesheets. None of 
the available literature has considered adhesively 
bonded honeycomb structures. Furthermore, the 
majority of the reported work has neglected either 
radiation or gas conduction heat transfer. Swann and 
Pittman used a finite difference model to study 
combined conduction and radiation in honeycomb 
structures and used then results to derive a semi- 


empirical relationship for the effective thermal 
conductivity as a function of geometric parameters and 
material properties. Stroud measured effective diermal 
conductivities of four honeycomb core panels over a 
temperature range of 670 to 1050 K and showed a root 
mean square deviation of 7% between his measured 
values and predictions using die Swann and Pittman 
semi -empirical relationship. Eftekliar, et al., 3 
developed a simplified analytical solution for two- 
dimensional conduction heat transfer hi a rectangular 
enclosure, neglecting radiation heat transfer. 

Fairbanks 4 formulated effective thermal conductivity of 
square cells for heat flow lateral to die principal axes of 
die cells for honeycomb core mirrors in space-based 
applications. Jones 5 modeled combined solid 
conduction and radiation heat transfer dirough 
evacuated honeycomb cored panels. Copenhaver, et 
al., 6 used die finite element numerical technique to 
model combined conduction/radiation in honeycomb 
core panels and used parameter estimation techniques 
in conjunction witii experimental heating data to 
estimate specific heat of the facesheets, emissivity hi 
die core and die conduction area of die core. 

Most of the reported work has concentrated on 
conduction and radiation as die prevalent modes of heat 
transfer and has neglected natural convection. 

Edwards, et al., 7 studied natural convection in 
honeycomb structures and determined die critical 
Rayleigh number for die onset of natural convection hi 
square cell honeycombs. The maximum calculated 
Rayleigh number for the titanium honeycomb core 
panel studied here was determined to be 7900, which 
was lower than the critical Rayleigh number for diis 
geometry, 61000, and therefore, natural convection was 
not considered as a mode of heat transfer in the present 
study. 

Honeycomb Core Panel Description 

The honeycomb core panel used hi this study was 
made of titanium 3A1-2.5V widi a density of 94.5 
kg/m 3 . The core was 25.4 mm thick and consisted of 
4.76 mm non-corrugated hexagonal cells widi 0.035 
mm foil thickness. The facesheets were made of 
titanium 6-22-22 and were 1.6 mm thick. The 
honeycomb core was attached to the facesheets ushig a 
modified epoxy adhesive film (FM5) supported by a 
polyester knit fabric. Two layers of the adhesive at 
0.49 kg/m 2 were applied to each facesheet. The 
adhesive thickness was determined to be 0.94 mm. 
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Swann and Pittman Model 


Swann and Pittman developed a semi-empirical 
model for heat transfer through honeycomb sandwich 
panels. 1 They assumed that the facesheets' thicknesses 
were very low to make the thermal contribution of the 
facesheets negligible, and thus simply modeled heat 
transfer through the honeycomb core structure. They 
considered solid conduction through honeycomb cell 
walls, gas conduction in the honeycomb cells, and 
radiation heat transfer hr tire honeycomb cell enclosure. 
They developed an empirical relationship for modeling 
radiation in the honeycomb core enclosure, and used a 
parallel thermal network model for modeling solid and 
gas conduction through die honeycomb core. The 
effective thermal conductivity, lq,, is given by: 


. . AA A/V 

k e = k f + k g ( 1 - — — ) + k,. 


( 1 ) 


where kf is the tiiermal conductivity of honeycomb core 
foil parent material, k g is the diermal conductivity of the 
gas contained in die panel, and AA/A is die ratio of die 
cross sectional areas of the solid core to die overall 
honeycomb cell. The radiation effective conductivity, 
k r , is given by: 

k r = 0.664(X + o.3) -0 ' 69 e L63(X - +1) "° S9 

L , <t(T l +T 2 )(T x 2 +T 2 2 ) 


where ct is the Stefan-Boltzmann radiation constant, 
and Tf and T 2 are the temperatures of the facesheets, L' 
is the honeycomb core height, X is the ratio of height to 
diameter for a cell in the honeycomb panel, and s is 
emissivity. It is assumed that the facesheets and 
honeycomb core have the same emissivity. The 
radiation effective conductivity function is an empirical 
relationship developed by Swann and Pittman using a 
trial and error process by comparing results with a finite 
difference model of the combined conduction/radiation 
heat transfer process in the honeycomb cell. 


equivalent diameter was therefore found to be 9.096 
mm. This would result in die height to diameter ratios, 
X, of 2.79 when ignoring die adhesive layers, and 2.58 
when including die adhesive layers. 


Experimental Measurements 

The effective thermal conductivity of the 
honeycomb core panel had been measured using drree 
different techniques at four diermophysical property 
measurement laboratories. Two laboratories utilized 
die guarded hot plate technique, 8 which is an absolute 
steady state technique. This technique is reported to 
have uncertainties of ±5% when measuring diermal 
conductivity of high- density rigid specimen. 8 One 
laboratory used the heat flow meter apparatus based on 
American Society of Testing and Materials (ASTM) 
Standard C5 1 8, 9 which is a steady-state comparative 
technique, requiring specimens of known transmission 
properties to calibrate die heat flux transducers used hi 
die apparatus. The measurement laboratory stated die 
general overall uncertainty of dieir measurement as 
±5%. One laboratory used a transient radiant step 
heating technique. One face of die specimen is 
subjected to a step radiant heating, and temperatures are 
measured on both sides of the specimen. The diermal 
diffusivity of the specimen is determined by 
minimizing the sum of square of differences between 
measured and predicted temperatures using a one 
dimensional finite difference model of the heat 
conduction equation through the sample assuming die 
test specimen is a homogenous material. The tiiermal 
conductivity is then calculated from the thermal 
diffusivity by knowing or measuring the specimen 
specific heat. No measurement uncertainties were 
available for the data generated using tiiis technique. 

All the measurements were performed on 203.2 mm 
diameter specimen, except for the radiant step heating 
technique which used a 203.2 mm diameter specimen 
for the room temperature measurement and a 50.8 mm 
diameter specimen for measurements above room 
temperature. 


The ratio of the cross sectional areas of the solid 
cores to the overall honeycomb cell required for 
calculations in Eq. (1) was calculated based on the ratio 
of the mass density of the honeycomb core to the bulk 
material, and for the current honeycomb geometry this 
ratio was 0.0208. For the Swann and Pittman semi- 
empirical relationship and the finite volume numerical 
heat transfer model, the hexagonal honeycomb cells 
were modeled as circular cells. The equivalent 
diameter of the circular cell was obtained by equating 
the perimeters of the hexagonal and circular cross 
sections as suggested by Swann and Pittman. The 


The variation of the measured thermal conductivity 
with average specimen temperature for the temperature 
range of 300 to 450 K using the measurements from die 
titt'ee different techniques at the four laboratories are 
shown in Figure 1 . The prediction of the effective 
diermal conductivity from the semi -empirical Swann 
and Pittman 1 model is also included in the figure. Error 
bars used for the steady-state techniques correspond to 
die stated ±5% measurement uncertainty; no error 
estimates were available for the radiant step heating 
technique. The first observation is that the 
measurements using die three techniques were 
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significantly different from each other. Even the 
measurements using the guar ded hot plate technique, 
which is an absolute measur ement technique and is the 
most frequently used method for thermal conductivity 
measurement, at the two different laboratories were 
significantly different from each other. The guarded 
hot plate data from the two laboratories agreed at room 
temperature, but diverged from each other as tire 
specimen temperature increased. Tire heat flow meter 
data were lower than the guarded plate data, but did not 
exhibit tire same trend as the other data. Tire thermal 
conductivity increased with increasing temperature up 
to 370K, and then asymptotically approached a constant 
value. The measur ements at this laboratory used tire 
National Institute of Standards and Technology (NIST) 
standard reference material (SRM) 1450b, a fibrous 
glass board, for calibrating tire heat flux transducers. 
Since the standard reference material had a lower 
thermal conductivity than tire honeycomb core parrel, 
the results were obtained by extrapolating tire heat flux 
transducers’ calibration curves which may have resulted 
in higher uncertainties than tire stated general 
uncertainty of ±5%, especially at higher temperatures. 
The radiant step heating technique produced tire lowest 
measured thermal conductivities and was closest to the 
Swann and Pittman predictions. 

Tire second observation is that the data for tire four 
measurements were all significantly different from the 
Swann and Pittman predictions. The average percent 
difference over the entire temperature range, with the 
corresponding 95% confidence intervals, 10 between the 
measurements and the Swann and Pittman predictions 
were 51.713.1, 60.9128.7, 28.4115.8, and 17.4128.4 for 
the two guarded hot plate data, heat flow meter data, 
and the radiant step heating data, respectively. The data 
are also presented in Table 1. 


measurement 

technique 

average percent 
difference with respect 
to Swann & Pittman 

Guarded hot Plate- 1 

51.7 ±3.1 

Guarded hot Plate-2 

60.9 ±28.7 

heat flow meter 

28.4 ±15.8 

radiant step heating 

17.4 ±28.4 


Table 1. Average percent difference of various 
experimental data with respect to Swann and Pittman 
model. 1 


Effective 

thermal 

conductivity, 

W/(m.K) 



Temperature, K) 


Figure 1 . Comparison of the measured effective 
thermal conductivities of the adhesively bonded 
honeycomb core panel using three techniques at four 
laboratories with Swann and Pittman model. 1 


Combined Conduction/Radiation Heat Transfer in 
Honeycomb Panel 

A steady state finite volume numerical method was 
used to solve the combined conduction/radiation heat 
transfer in the honeycomb core panel. The governing 
equations were applied to a volume representing one 
cell of the honeycomb core and the associated 
facesheets and adhesive layers. The hexagonal core of 
the honeycomb was approximated as having a circular 
cross section. Swann and Pittman had used the same 
approximation in their finite difference numerical 
solution of the problem. This approximation reduced 
the problem to an axisymmetric case, hi the absence of 
natural convection, the governing partial differential 
equation for conduction heat transfer is: 


1 d , , dT. 9 ,, dT, 

r dr dr dz dz 


The boundary conditions are: 


k— (L,r) = q' 
dz 

T(0, r) =T, 

— (z,0) = 0 
dr 

Tt 1 

— (z,ii) = 0 
dr 


(3) 


(4a) 

(4b) 

(4c) 

(4d) 


The first two boundary conditions represent the 
boundary conditions applied to the two facesheets of 
the honeycomb core panel. The first boundary condition 
represents an applied constant heat flux, while the 
second boundary condition represents a constant 
temperature boundary condition. The last two 
boundary conditions represent symmetry conditions at 
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the cell center and boundary. A schematic of the 
geometry considered for this simulation is shown in 
Figure 2. The finite volume formulation of the 
conservation of energy equation was developed by 
performing an energy balance at each node: 11 

Uy(T if j-Ti_i J ) + v y (T y -Ti + i J ) + 

W y ( T iJ - T iJ-i > + X i,j ~ T i.i + i > - ( 5 ) 

E i,j+ G i,j=° (i = l,nj;j = l,nj) 

where Ty was tlie nodal temperature, Uy, Vy, Wy, and 
Xy were the thermal conductances, which in their most 
general form for non-uniform nodal spacing (nodes 
being at the centroids of the volume cells) and for nan- 
homogeneous material were given by: 


u y 

_ kyArj +k i _ 1J Ai-_ 1 

(fi -0.5Ari)Azj 

(6a) 

Ah +Ar 1 _ 1 


v « 

_ kyAq + k i+l j Ar i+1 

(r ; +0.5Ar i )Azj 

(6b) 

Ati+Ar i+ i 

r i+i - r i 

w i,. 

kt.jAzj +k iJ _ 1 Az J _ 

■i f i Ar i 

(6c) 

’ Azj+Azj.j 

z r z H 

Xy 

ki.jAzj +kyyAz j+l ,:Ar 

(6d) 

AZ ; + AZ : , , 

z 14.1 - Z i 


where r ; and Zj were the radial and axial nodal positions, 
and Ay and Azj were the dimensions of the control 
volume cells in the radial and axial directions, 
respectively, n; was the total number of cells in the 
radial direction, while iij was the total number of cells 
in the axial direction. The first ratio term on the right 
hand side of Eq. (6a) through (6d) used the rule of 
mixtures to find the effective local thermal conductivity 
between two adjacent cells composed of different 
materials, and reduced to the thermal conductivity of 
the material if two adjacent cells were made of the same 
material. The numerator of the second ratio term on the 
right hand side of Eq. (6a) through (6d) was the 
effective area of the finite volume cell face for 
conduction into the cell, while the denominator was the 
distance between two adjacent cell centroids. The 
equations were modified for the cells on the boundaries 
based on the corresponding boundary conditions. The 
term By in Eq. (5 ) represented the external heating 
applied at one facesheet, therefore it was zero 
everywhere in the domain except at the boundary nodes 
at z=L (i=l,iti, j=nj), and was obtained by multiplying 
the applied heat flux by the cross sectional area of the 
respective control volume cells. The term Gy in Eq. (5) 
represented radiant heat in the honeycomb cell 
enclosure and was zero everywhere except for control 


volume cells with control surfaces that participated in 
radiation exchange in the enclosure. These consisted of 
the control surfaces on the honeycomb core and on the 
adhesive layers at the top and bottom of honeycomb 
cores. Therefore, the term Gy was calculated by 
multiplying the radiant heat flux, q", by the 
corresponding control volume surface areas. 


A 





facesheet 
p adhesive 

|“ 


-foil 


r r, 


Figure 2. Schematic of honeycomb core geometry. 


The radiant heat fluxes were calculated at each 
iteration by performing radiation exchange analysis in 
the enclosure assuming gray surfaces, with each control 
volume cell assumed to be isothermal with uniform 
surface radiosity, and furthermore, that both emitted 
and reflected radiation were diffuse. Letting 0 k 
represent nodal temperature Ty mapped to a single 
index, the incident radiant heat flux at each control 
volume surface in the honeycomb enclosure involved in 
radiation exchange was calculated from: 12 


where: 


q k ,i = Z A ik o0 k , 1 < i < N 

k=l 

(7) 

A* =-^(8*-**) 
l-6i 

(8a) 

Ta 

e 

ii 

rS 

(8b) 

^ 5 ik -(l-6 i )F A A k 

®ik= — 

(8c) 


where N was the total number of control volume 
surfaces involved in radiation exchange hi the 
enclosure, Bi was the surface emissivity for each control 
volume cell. 5 lk was the Kronecker delta, and FAi-Ak 
was the shape factor for radiation exchange between 
surfaces designated as Aj and A k . The shape factors for 
the radiation exchange were calculated using the shape 
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factor equations provided by Swann and Pittman 1 for 
die same geometry. 

The modified Newton-Raphson method was used 
to solve the nonlinear set of equations: 13 

[J] m {A0} m+1 =-{F}"' (9a) 

{0} w+1 ={6}" +{A0} W+1 (9b) 

where die superscript m denotes die m' h iteration, [J] m 
was die Jacobian matrix, { A0 } m+1 was die nodal 
temperature increments vector caused by the residual 
load vector (F} m . The Jacobian matrix was defined as: 


J 


i.j 


3 ^ 

30; 


( 10 ) 


The Jacobian matrix was calculated by finite-difference 
numerical differentiation at die first iteration and was 
held constant. Once die numerical solution had 
converged, die effective thermal conductivity was 
calculated based on die applied heat flux boundary 
condition and die converged average temperature 
difference between die two facesheets and the sample 
overall diickness using die Fourier’s law of heat 
conduction. 14 


Published results were used for die thermal 
conductivity of the titanium honeycomb core (Ti-3A1- 
2.5 V), 15 and air. 14 No data were found for the 
facesheets material, Ti-6-22-22. Therefore thermal 
properties for Ti-6A1-4V were used for the facesheets. 16 
The thermal conductivity of the adhesive layer was 
unknown. It was decided to investigate two values for 
die thermal conductivity of the adhesive that would 
provide an upper and lower range for this property. It 
was assumed that the ratio of the thermal conductivity 
of the adhesive to the thermal conductivity of air over 
the temperature range of interest was either 10 or 100. 
This would result in thermal conductivity of the 
adhesive ranging between 0.25 1 and 2.5 1 W/niK at 
room temperature. This would cover the range of 
thermal conductivity values for typical adhesives and 
epoxies. The emissivity of titanium foil and adhesive 
layer were assumed to be independent of temperature in 
the temperature range of interest and equal to 0.3 17 and 
0.8, respectively. The emissivity of nomnetallic solids 
generally varies between 0.7 and 0.9, 12 therefore, an 
emissivity of 0.8 was assumed for the adhesive layer. 


Modified Swann and Pittman Model 


It was further decided to investigate a modification 
to the original Swan and Pittman semi-empirical 
relationship. The proposed modification incorporated 
the facesheets and the adhesive in the thermal model. 
The equivalent thermal network model consisted of five 
thermal resistances in series arrangement. Two of the 
resistances were for the solid conduction through the 
facesheets. Two other resistances were the effective 
resistances for the parallel arrangement of solid 
conduction through the adhesive and the segment of 
honeycomb foil covered by the adhesive. The last 
resistance was the effective resistance for the parallel 
arrangement of conduction through the gas (air) and 
remaining honeycomb foil along with the superimposed 
radiation conduction term. Titus, the effective thermal 
conductivity for this thermal network was obtained 
from the following relationship: 



-'ad 


AA , AA 

-r+ k ad(i— — ) 
A A 


+ 


V 

AA AA 

f X +kg( 1 ~X )+kl ' 


(ID 


where Lf s , L ad , and L were the thicknesses of the 
facesheet, adhesive, and the honeycomb core panel, 
respectively, V was the diickness of the gas column in 
the honeycomb core (L'=L-2L fs -2L ild ), and kf s , kf, k ad , 
and k g were the thermal conductivities of the facesheet, 
foil, adhesive and gas, correspondingly. Note that k r 
was the radiant thermal conductivity calculated from 
Eq. (2). The original Swann and Pittman model was 
developed assuming that the facesheets and foils had 
the same emissivity. The adhesive layers at the top and 
bottom of the honeycomb cell had a different emissivity 
value than the honeycomb core foil. In order to utilize 
the Swann and Pittman radiation model, a weighted 
emittance based on the ratio of the surface area of each 
material to the total surface area involved in the 
radiation exchange in the enclosure was calculated. 

This yielded a weighted emissivity of 0.382. 


Results and Discussion 


Before proceeding to presentation of the finite 
volume numerical results for the adhesively bonded 
geometry, the finite volume numerical model was 
applied to the honeycomb core panel geometry without 
any adhesive layers in order to assess the accuracy of 
the numerical results by comparison with the Swann 
and Pittman predictions using Eq. (1). It was assumed 
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ill the numerical model that the emittance of the 
facesheets was equal to the emittance of the honeycomb 
foils, 0.3. The effective thermal conductivities of the 
honeycomb core panel without the adhesive layers at 
various temperatures were calculated using the 
converged finite volume numerical solution of the 
combined conduction and radiation heat transfer. The 
results and die Swann and Pittman predictions are 
shown in Figure 3, where the effective thermal 
conductivities are plotted versus average honeycomb 
panel temperatures. The Swann and Pittman model 
overpredicted die numerical results below 300 K and 
underpredicted die numerical results above 300 K, widi 
die difference increasing widi increasing temperature. 
The average difference between die Swann and Pittman 
model and die finite volume numerical predictions was 
— 3.5±6.6 percent. Considering diat die Swann and 
Pittman model ignores die facesheets and provides a 
semi-empirical relationship for modeling die radiation 
heat transfer, diis close agreement between die finite 
volume numerical predictions and die semi-empirical 
results validated die numerical mediod. These results 
were also in general agreement widi results from 
Stroud ’ who had showed a root mean square deviation 
of 7% between his experimental values and predictions 
using die Swann and Pidman semi-empirical 
relationship over a temperature range of 670 to 1050 K. 


Effective 

thermal 

conductivity, 

W/(m.K) 



Temperature, K 


Figure 3. Comparison of finite volume numerical model 
and Swann and Pittman model 1 for honeycomb core 
panel without adhesive. 


conductivity values is presented in Figure 4. The 
experimental data and die Swann and Pittman 
predictions are also shown in die figure. The average 
percent differences for die effective thermal 
conductivity of the Swann and Pittman model and die 
four experimental sets of data widi respect to die 
average numerical predictions are presented in Table 2. 
The average numerical prediction was obtained by 
averaging die data obtained at die two different values 
of adhesive thermal conductivity. It is evident diat die 
radiant step heating data had provided die best 
agreement widi die numerical predictions, widi die data 
clearly falling in between die two curves corresponding 
to die two diermal conductivity values of the adhesive. 
Using die average numerical predictions as die basis for 
reference, die percent difference for die radiant step 
heating measurements was 1.7±5.9; the percent 
difference for die heat flow meter measurements was 
12.9+26.8; while die two sets of guarded hot plate data 
produced differences larger tiian 38 percent. These 
results do not in any way question die validity of the 
steady state techniques, but imply diat die 
measurements obtained for die samples were 
inaccurate. 


Effective 

thermal 

conductivity, 

W/(m.K) 



Temperature, K 


Figure 4. Comparison of effective thermal conductivity 
predicted using the numerical model with previous 
measurements and the Swann and Pittman model 1 for 
thermal conductivity of adhesive to ah' ratios of 10 and 


100 . 


The effective thermal conductivities of the 
honeycomb core panel with the adhesive layers at 
various temperatures were calculated using the 
converged finite volume numerical solution of the 
combined conduction and radiation heat transfer. As 
discussed previously, results were obtained with the 
ratio of the thermal conductivity of adhesive to ah of 10 
and 100. The variation of effective thermal 
conductivity with temperature calculated horn the 
numerical results at the two adhesive thermal 
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measurement 
technique or 
prediction method 

average percent 
error 

Guarded Hot Plate- 1 

38.0 ± 18.1 

GHP-2 

47.3 ±11.1 

heat flow meter 

12.9 ± 26.8 

radiant step heating 

1.7 ±5.9 

Swann & Pittman 

-15.2 ±12.6 


conductivity of adhesive to air ratio of 100 was 
2.3±1 1.6. These percent differences were significantly 
lower compared to the Swann and Pittman model and to 
all the steady state experimental results. It can be 
concluded that the modified Swann and Pittman model 
is a good first-order approximation for predicting die 
effective thermal conductivity of adhesively bonded 
honeycomb core panels. 


Table 2. Average percent differnce of various 
experimental data and Swann and Pittman model 1 widi 
respect to average predictions by the finite volume 
numerical model. 


The percent difference between die average 
numerical results and die Swann and Pittman model 
was — 15.2±12.6. It is obvious diat die Swann and 
Pittman semi-empirical relationship is insufficient for 
modeling die heat transfer in die adhesively bonded 
honeycomb panel used in this investigation. It can be 
concluded diat die adhesive layer has a significant 
effect on die overall thermal performance of 
honeycomb core panels and should not be ignored in 
die calculations, hi a honeycomb core panel at 
moderate temperatures where radiation is not a 
significant mode of heat transfer, most of die heat 
transfer takes place dirough the honeycomb core foils. 

In the absence of the adhesive layers, heat is conducted 
from the facesheets to the foil through the small cross 
sectional areas of the foil, JtDtf, where tf is foil 
diickness and D the hydraulic diameter. In the presence 
of adhesive layers heat is also transferred front the 
facesheets to the adhesive layer and then front the 
adhesive layer to the foil through an area equal to 
7tDL a( j, where L ad is the adhesive thickness. Depending 
on die diermal conductivity and thickness of the 
adhesive layer, this added path for conduction of heat 
from the facesheet to foil through the adhesive, could 
result in higher heat conduction rates through the 
overall honeycomb core panel. 

For the honeycomb core panel geometry of the 
present study, the finite volume numerical results and 
the modified Swann and Pittman results for the thermal 
conductivity of adhesive to air ratio of 10 are shown in 
Figure 5. The Swann and Pittman predictions are also 
included in the figure. The modified Swann and 
Pittman overpredicted the numerical results below 400 
K and underpredicted the numerical results above this 
temperature, but it is consistently higher compared to 
the original Swann and Pittman. The percent difference 
between the modified model and numerical results was 
3.2±12.5. The percent difference between die modified 
model and finite volume numerical results for thermal 


Effective 

thermal 

conductivity, 

W/(m.K) 



Temperature, K 


Figure 5. Compar ison of effective thermal conductivity 
predicted using die modified Swann and Pittman model 
with numerical model and die Swann and Pittman 
model for diermal conductivity of adhesive to air ratios 
of 10. 


The finite volume numerical model was used to 
paramed ical ly study the effects of adhesive thickness 
and thermal conductivity on the effective thermal 
conductivity at a mean temperature of 300K. The ratio 
of the adhesive layer thickness (on each side) to 
honeycomb foil height was varied between zero and 
0.05, with zero representing a construction widi no 
adhesive. The geometry studied in diis investigation 
had an adhesive thickness to foil height ratio of 0.037. 
Data were obtained for the ratio of thermal conductivity 
of adhesive to air being equal to 1, 5, 10, 25, 50, and 
100. All the other parameters were kept fixed at die 
values corresponding to the geometry studied hi diis 
investigation. Results presented in Figure 6 illustrate 
the variation in effective thermal conductivity increase 
due to the adhesive with the ratio of the adhesive layer 
diickness to honeycomb foil height for various values 
of the ratio of the adhesive to air thermal conductivity. 
All the data converged to zero percentage increase 
when the adhesive thickness was zero. The effective 
thermal conductivity increased with increasing adhesive 
diickness and adhesive thermal conductivity. 
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Percent 
increase in 
effective 
thermal 
conductivity 



Ratio of adhesive thickness to foil height 


Figure 6. Percent increase in effective thermal 
conductivity as a function of ratio of adhesive thickness 
to foil height for thermal conductivity of adhesive to ah 
ratios of 1, 5, 10, 25, 50, and 100. 


neglecting the influence of the adhesive in the overall 
heat transfer. 


Concluding Remarks 

The effective thermal conductivity of an adhesively 
bonded honeycomb core structure had been measured 
using three different techniques at four thermophysical 
property measurement laboratories over the temperature 
range of 300-500 K. Two laboratories had utilized the 
guarded hot plate technique, one laboratory had used 
the heat flow meter apparatus based on ASTM Standard 
C518, which is a steady-state comparative technique, 
and one laboratory had used the transient radiant step 
heating technique. The four sets measurements varied 
significantly from each other and from the semi- 
empirical Swann and Pittman model. 1 


The data for adhesive to air thermal conductivity 
ratio of unity illustrated die influence of adhesive 
emissivity and thickness on die radiation heat transfer. 
The adhesive did not change die conduction heat 
transfer through die geometry since it had die same 
diermal conductivity as air, and dius die effect of the 
adhesive on die radiation heat transfer could be 
isolated. The adhesive had a higher emissivity than die 
facesheets, and the adhesive thickness changed the 
overall height of the enclosure for radiation exchange. 
The effective conductivity increased with increasing 
adhesive thickness and resulted in a maximum three 
percent higher effective thermal conductivity at the 
adhesive thickness to foil height ratio of 0.05. 

The data for adhesive thermal conductivities larger 
dian that of air exhibited a larger increase in the 
effective thermal conductivity with increasing adhesive 
diickness and with increasing adhesive thermal 
conductivity. The higher increases were due to the 
combination of increases in the radiation heat transfer 
and solid conduction heat transfer. As discussed 
previously, the adhesive layer provided a larger surface 
area with higher thermal conductivity to conduct heat 
from the facesheets to the foils than did the foil edge. 
The difference between the data for adhesive thermal 
conductivity ratios greater than one compared to ratios 
equal to one provided a rough estimate of the 
magnitude of the higher heat transfer due to conduction 
within the adhesive. At an adhesive to air thermal 
conductivity ratio of 100, the adhesive’s thermal 
conductivity approached that of the titanium facesheets, 
and the problem reduced to heat transfer in a 
honeycomb structure with an effectively thicker 
facesheet and smaller foil height. The data in the figure 
provide an estimate of the errors encountered when 


The combined steady-state radiation and 
conduction heat transfer in the honeycomb core panel 
was modeled using a finite volume numerical 
formulation. The nonlinear governing equations were 
solved using a modified Newton-Raphson method. The 
effective thermal conductivities were calculated from 
the converged steady-state solutions of the numerical 
finite volume model using thermal conductivity of 
adhesive to air ratios of 10 and 100. The transient 
radiant step heating experimental data were found to be 
in best agreement with the numerical finite volume 
predictions with a percent difference of 1 .7±5.9. 

A modification of the Swann and Pittman semi- 
empirical relationship which incorporated the 
facesheets and adhesive layers in the thermal model 
provided satisfactory agreement with the finite volume 
numerical results. Tire percent difference between the 
numerical predictions and the modified model for 
thermal conductivity of adhesive to ah' ratios of 10 and 
100 were 3.2±12.5 and 2.3±1 1.6, respectively. It can 
be concluded that the modified Swann and Pittman 
model is a good first-order approximation for predicting 
the effective thermal conductivity of adhesively bonded 
honeycomb core panels. A parametric study was 
conducted to determine the effects of adhesive 
thickness and thermal conductivity on the effective 
thermal conductivity. It was determined that the 
effective thermal conductivity of the panel increases 
with increasing adhesive thickness and adhesive 
thermal conductivity. 
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